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ABSTRACT 

We examine the energy-resolved pulse profiles of 51 flux oscillations observed during the decline of 
thermonuclear X-ray bursts from accreting weakly-magnetized neutron stars with the Rossi X-ray Timing 
Explorer. We find that the fractional rms amplitudes of the oscillations increase as a function of energy 
by 0.25% keV" 1 to 0.9% keV" 1 between 3-20 keV, and are as large as 20% in the 13 - 20 keV band. We 
also show that the pulses observed in the higher energy bands generally lag behind those in lower energy 
bands by 0.002 cycles keV -1 to 0.007 cycles keV -1 between 3 — 20 keV. This amounts to total delays of 
0.03-0.12 cycles between the lowest and highest energy bands, or time delays that range from 100-200 
/is. We then model the oscillations as flux variations arising from temperature patterns on the surfaces 
of rapidly rotating neutron stars. In this framework, we find that the increase in the pulse amplitude 
with photon energy can be explained if the cooler regions on the neutron star emit in the lower energy 
bands, reducing the flux variations there. On the other hand, the Doppler shifts caused by the rapid 
rotation of the neutron star should cause the hard pulses to precede the soft pulses by about 0.05 cycles 
(100 /is), in contrast to the observations. This suggests that the photons originating from the stellar 
surface are reprocessed by a hot corona of electrons before they reach the observer. 
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1. INTRODUCTION 

Thermonuclear bursts from accreting weakly- 
magnetized neutron stars provide an excellent opportu- 
nity to study emission originating from the stellar surface 
(see Lewin, van Paradijs, & Taam 1993 for a review). 
For instance, the unstable helium burning that produces 
the bursts may leave asymmetries in the surface bright- 
ness of the star that can be observed as pulsations at 
the stellar spin frequency (e.g., Fryxell & Woosley 1982; 
Strohmayer et al. 1997). Indeed, oscillations that are ap- 
parently produced by bright regions on rapidly rotating 
neutron stars now have been observed during bursts from 
ten neutron star low-mass X-ray binaries (LMXBs) with 
the Rossi X-ray Timing Explorer (RXTE; see Strohmayer 
2001, for a review). Their frequencies range between 270- 
620 Hz (Muno et al. 2002a), which are consistent with 
the assumption that these LMXBs are the progenitors of 
recycled millisecond radio pulsars (Alpar et al. 1982; Rad- 
hakrishnan & Srinivasan 1982). Most importantly, after a 
small frequency drift is accounted for, the oscillations are 
extremely coherent. The oscillations have quality factors 
Q = v / Av > 1000 during individual bursts (Strohmayer 
et al. 1996; Strohmayer k, Markwardt 1999, 2002; but see 
Muno et al. 2002a), and they appear at frequencies that 
are stable to a few parts in 1000 in bursts observed over 
several years from any given source (Strohmayer et al. 
1998b; Giles et al. 2002; Muno et al. 2002a). The stability 
and coherence of these oscillations provides compelling 
evidence that they originate from brightness patterns on 
rotating neutron stars. 

Burst oscillations therefore allow one to study a range 
of questions, including the physics of nuclear burning on 
the stellar surface, the effects of gravitational self-lensing 
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by the neutron star, and the equation of state of dense 
nuclear matter. The most constraining data so far have 
come from the amplitudes and the profiles of the oscilla- 
tions. The fractional rms amplitudes of the oscillations 
can be as large as 50% during the first second of a burst 
(Strohmayer et al. 1998a). If the oscillations result from 
two antipodal bright regions on the surface, as has been 
proposed based upon both observational and theoretical 
considerations (Miller, Lamb, & Psaltis 1998; Miller 1999), 
then amplitudes this large can only be produced if the 
neutron star is larger than R/Rg c h = 3.1 (Miller & Lamb 
1998; Nath, Strohmayer, & Swank 2002), where R Sc h is 
the Schwarzschild radius. Moreover, the profiles of the 
oscillations are nearly perfectly sinusoidal — the ampli- 
tudes of the harmonic components are less than 0.5% of 
those of the main signals. This suggests that the patterns 
forming on the stellar surfaces are either highly symmetric 
or restricted to near the rotational poles (Muno, Ozel, & 
Chakrabarty 2002b). 

The surface velocity of the neutron star (ss 0.1c for a 
10 km star spinning at 500 Hz) can be measured from 
variations in the amplitude and phase of the oscillations 
as a function of energy, which would provide another con- 
straint on the stellar radius. The amplitudes of the os- 
cillations should increase as a function of energy because 
Doppler shifts cause the Wien tail of the spectral distribu- 
tion to move in and out of the highest energy bands (Miller 
& Lamb 1998). However, the amplitudes will also be set 
by the surface temperature distribution, since a finite sur- 
face temperature produces a constant background against 
which a hot region is less distinct at low energies (Page 
1995). The relative phases of the waveforms as a func- 
tion of energy, on the other hand, are affected only by the 
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Fig. 1. — The energy dependence of the amplitudes of the averaged oscillation profiles from six sources. For each source, oscillations were 
only averaged for a single epoch during which the gain of the PCA was relatively constant. In all cases, the amplitudes increase significantly 
as a function of the photon energy. 



Dopplcr shifts of the tail of the spectrum, which move the 
peak of the high-energy pulse to an earlier phase (Braje, 
Romani, & Rauche 2000; note that Doppler beaming and 
photon arrival time delays also shift the overall phase of 
the entire pulse). Evidence that the hard pulse leads the 
soft has been reported in a burst oscillation from Aql X-1 
by Ford (1999), although further studies question whether 
this is due to rotational motion (Fox 2000). 

In this paper, we present a study of the amplitudes and 
phases of burst oscillations as a function of photon energy. 
In Section 2 we analyze the energy-resolved waveforms of 
59 oscillations observed with the Proportional Counter Ar- 
ray (PCA) on RXTE. In Section 3, we present theoreti- 
cal waveforms produced from the models of Muno et al. 
(2002b), taking into account the response of the PCA to 
facilitate comparison with the data. We consider several 
simple temperature distributions on the neutron star. In 
Section 4, we explore whether we can measure Doppler 
effects in the observed waveforms. 



2. OBSERVED WAVEFORMS 

As part of an ongoing study of thermonuclear bursts ob- 
served with the RXTE Proportional Counter Array (PCA; 
Jahoda ct al. 1996), we identified 68 X-ray bursts contain- 
ing millisecond oscillations in data from eight sources that 
were in the public domain as of September 2001 (Muno et 
al. 2002a). The PCA is composed of five gas-filled propor- 
tional counter units (PCUs) that are sensitive to X-rays 
between 2.5-60 keV with an energy resolution of about 
15%. Between 2 and 5 PCUs were active during the ob- 
servations in our sample. Each PCU has an effective area 
of 1200 cm 2 and is capable of recording photon arrival 
times with microsecond time resolution. For most of the 
observations, data with both 2~ 13 s (122 /zs) time reso- 
lution and between 8-64 energy channels were recorded 
on an event-by-event basis, and were used to study the 
profiles of the oscillations. 

We first modeled the frequency evolution of the oscilla- 
tions observed in the entire PCA bandpass using a phase 
connection technique similar to that used in pulsar tim- 
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ing (Manchester & Taylor 1977; Muno et al. 2002a). To 
implement this technique, we folded the data in short in- 
tervals (0.25-0.5 s) using a trial phase model, which was 
then refined through a least x 2 nt to the residuals. This 
provides excellent frequency resolution and a statistical 
measure of how well the model reproduces the data. We 
have previously used the resulting phase models to study 
the frequency evolution, amplitudes, and profiles of the 
burst oscillations (Muno et al. 2002a; Muno et al. 2002b). 

We then applied the best-fit phase models to produce 
folded profiles in 7 energy channels. Since the gains of 
the proportional counters have been changed several times 
during the mission, the energies corresponding to the de- 
tector channels vary. Data with sufficient energy and time 
resolution were only used to record burst oscillations in 
gain epochs 3 (1996 April 15 through 1999 March 22), 
4 (1999 March 22 through 2000 May 13), and 5 (2000 
May 13 to present). We were unable to choose detector 
channels with identical energy boundaries in different gain 
epochs, since the high time resolution data were usually 
taken with too few energy channels. Therefore, we defined 
the energy boundaries by detector channels 5, 13, 18, 28, 
36, and 53. These boundaries roughly correspond to <2, 
5, 7, 10, 13, and >20 keV during gain epoch 3 (see also 
Table 1). The low (< 2 keV) and high (> 20 keV) energy 
channels generally were dominated by background events, 
and were therefore ignored during our analysis. We found 
that 51 oscillation trains from 6 sources were sampled with 
adequate time and energy resolution for this study. The 
sources, oscillation frequencies, and numbers of oscillations 
examined are listed in Table 2. 

We note that the energy-resolved data that we used (Ta- 
ble 2) contained up to 30% fewer counts than the pro- 
files created for a single energy channel in Muno et al. 
(2002b). This is because the RXTE data buffers could 
only record energy-resolved data for a fraction of each sec- 
ond before they were filled during the brightest segments of 
many bursts. Other modes designed specifically to record 
bursts were available for producing the frequency mod- 
els (so-called burst-catcher modes), but sacrificed energy 
resolution in favor of high time resolution. The burst- 
catcher modes were not used to produce energy-resolved 
waveforms. 

We measured the amplitudes of the oscillations as a 
function of energy by computing a Fourier power spectrum 
of the folded profiles in units of total counts per phase bin. 
We normalized the power according to Leahy et al. (1983), 
so that the fractional rms amplitude at any multiple of the 
oscillation frequency is 



X=[ T 



1/2 



(1) 



T — R 

where P n is the power at the nth bin of the Fourier spec- 
trum, I 7 is the total number of counts in the profile, and 
B 1 is the estimated number of background counts in the 
profile. The background was taken to be the persistent 
emission during the 15 s prior to the burst. This equa- 
tion is valid so long as the phase and frequency of the 
oscillation is known, as it is by design for our folded pro- 
files. Uncertainties and upper limits on the amplitudes are 
computed taking into account the distribution of powers 
from Poisson noise in the spectrum, using the algorithm 

4 http : //heasar c . gsf c . nasa . gov/docs/sof tware/ f tools/ 



in the Appendix of Vaughan et al. (1994). We measured 
the phases of the oscillations by linear least-squares fits 
of sinusoids to the folded profiles. The uncertainties on 
the phases were determined from the diagonal values of 
the covariance matrices that were derived as part of the 
least-squares algorithm (Press et al. 1992). 

We also summed the energy-resolved profiles for all of 
the oscillations that occurred during an individual gain 
epoch from each source (Table 2). Although the gain of 
the PCA drifted even during a single epoch, we confirmed 
that the drift was not large enough to affect the summed 
profiles. We did so by folding a set of simulated profiles 
(see Section 3) through different responses generated by 
the FTOOL 4 pcarsp that spanned gain epoch 3, measur- 
ing the amplitude and phase of each individual oscillation, 
and comparing it to those of the summed profile. The 
values were identical to within less than 1%. 

2.1. Amplitudes as a Function of Energy 

First, in order to determine whether the oscillation am- 
plitudes in individual bursts change significantly as a func- 
tion of photon energy, we calculated the \ 2 statistic under 
the assumption that the amplitudes are constant: 



E 



Ai A-i 

OA, 



(2) 



where Ai is the fractional rms amplitude in each energy 
band, Ai is the mean amplitude, and a Ai is the uncer- 
tainty derived according to the method of Vaughan et al. 
(1994). We find that 34 of 51 oscillations are inconsis- 
tent with a constant amplitude as a function of energy at 
the 90% confidence level. The remaining oscillations that 
appear to have constant amplitudes typically have fewer 
counts in the folded profile: the mean number of counts 
in a folded profile is 5.5 x 10 4 , while those with constant 
amplitudes have on average 3.5 x 10 4 counts. 

To quantify how the amplitude varies with energy in 
these oscillations, we fit them with a linear function using a 
least-squares technique. A linear fit is only a crude approx- 
imation to the functional form expected from a hot region 
on a rotating neutron star (see Section 3), but it is ade- 
quate for making simple quantitative comparisons between 
our observational and theoretical results, given the limited 
statistics in most of our measurements. We found that 28 
oscillations were consistent with a linear increase in am- 
plitude as a function of energy. The average of the best-fit 
slopes for these oscillations was 0.8% keV -1 . The smallest 
total detectable increase in amplitude was about 2% be- 
tween the 2 keV and 20 keV bands, while the increase was 
as large as 30% in several oscillations from 4U 1636—536 
and 4U 1702—429. Only 6 oscillations were inconsistent 
with both constant amplitudes and linear trends. Two of 
these (occurring on 1996 December 21 17:30:34 and 1999 
June 21 19:13:02 from 4U 1636—536, where these and all 
subsequent times are in Barycentric Dynamical Time) ex- 
hibit random variations in amplitude, while the remaining 
four (occurring on 2000 June 15 05:13:02 and 2000 August 
9 09:00:34 from 4U 1636-536, and on 1997 September 20 
10:08:27 and 1999 August 19 12:13:22 from 4U 1728-34) 
show deviations from a linear trend that are consistent 
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Fig. 2. — The energy dependence of the phases of the averaged oscillation profiles from six sources (compare Figure 1). Phase zero is defined 
by the profile in the total PCA bandpass. A positive value of the phase lag indicates that the profile in that energy range arrives later than 
the average profile. 



with the functional form expected from a hot region on a 
rotating neutron star (see Section 3). 

We then performed the same analysis for 9 summed 
pulse profiles from 6 sources. We list the results in Table 3, 
and we display the amplitude as a function of energy in 
Figure 1. The amplitudes from all of the summed profiles 
increase as a function of energy with > 99.9% confidence. 
In most cases, the increase is consistent with linear func- 
tions with slopes between 0.5-0.8% keV -1 . This translates 
to a total increase in amplitude from 5-16% between 2- 
20 keV. Only the oscillations from 4U 1636—536 in gain 
epoch 5 are inconsistent with a linear trend. These have 
a higher amplitude in the 2-5 keV band than would be 
expected from extrapolating the trend at higher energies 
(Fig. 1), but are still roughly consistent with the functional 
form expected from a hot spot. 

2.2. Phase as a Function of Energy 

We then examined the relative phases of the oscilla- 
tions as a function of energy. The reference phases were 



taken to be those of the oscillations in the total 2-60 keV 
bandpass. An increasing value of the phase lag indicates 
that the pulse seen at high energies arrives later than 
(lags) that at low energies. Of the individual oscilla- 
tion trains, we find that only 13 exhibit phases that vary 
as a function of energy at the 90% confidence level. Of 
these 13 oscillations, 6 appear to change in phase with a 
slope of -0.005(3) to 0.019(5) cycles keV" 1 (1999 April 
29 01:49:24 from 4U 1636-536; 1999 April 10 09:54:12 
from MXB 1659-298; 1999 February 22 04:54:01 from 
4U 1702-429; 1999 August 19 15:50:57 from 4U 1728-34: 
and 1996 July 14 04:23:45 and 1999 February 26 17:10:56 
from KS 1731—260). Seven oscillations are not well-fit 
by linear trends, but exhibit random variations in phase 
as a function of energy (2000 August 9 01:22:23 and 
09:00:34 from 4U 1636-536; 1997 July 26 09:09:11 from 
4U 1702-429; 1997 September 21 18:10:31, 1997 Septem- 
ber 27 11:17:45, 1999 January 31 21:56:41, and 1999 Au- 
gust 19 12:13:22 from 4U 1728-34). The weighted average 
of the phase lag for the 44 oscillations that can be ade- 
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quately fit with a linear trend is 0.0024(5) cycles keV -1 . 

In Figure 2 we plot the phase as a function of energy 
for 9 summed profiles from 6 sources, and in Table 4 we 
compile the results. Applying the x 2 statistic to the data, 
we find that the phase varies as a function of energy at 
the 90% confidence level in 5 of the 9 oscillations. In 4 of 
these 5 oscillations, the phases are consistent with a linear 
increase as a function of energy with a slope of between 
0.002(1) and 0.007(3) cycles keV" 1 . Thus, the pulses at 
20 keV lag behind those at 2 keV by 0.04-0.14 cycles, 
which translates to time delays of 100 ^s (4U 1702-429) 
to 200 /us (MXB 1659-298). The summed oscillation from 
4U 1728—34 in gain epoch 3 appears to exhibit random 
variations in phase as a function of energy, as is evident in 
Figure 2. 

We note that our sample includes the oscillation from 
Aql X-l in which Ford (1999) reported that the pulse at 
high energies precedes that at low energies in the last half 
of the oscillation train. We find that these soft lags are sig- 
nificant only in the interval that was searched, while ear- 
lier the oscillation exhibited hard lags (see also Fox 2000) . 
Over the entire burst, the soft and hard lags cancelled, 
yielding no net delay between hard and soft photons. We 
searched for similar changes in the energy dependence of 
the phases of other oscillations by dividing the oscillations 
into two intervals containing an equal number of counts. 
We did not find any other significant examples of phase 
lags switching from positive to negative. 

3. MODEL PROFILES 

In this section, we examine how the rotation of the neu- 
tron star affects the energy dependence of the burst oscilla- 
tions. To that end, we produce theoretical energy-resolved 
light curves from circular hot spots on a rotating neutron 
star. We use the techniques outlined by Pechenick, Fta- 
clas, & Cohen (1983), Braje et al. (2000), and Weinberg, 
Miller, & Lamb (2001). The amplitudes and profiles of 
our theoretical oscillations in the entire PCA bandpass (2- 
60 keV) are discussed extensively in Muno et al. (2002b). 

The energy spectra observed during bursts can ade- 
quately be modeled as blackbody emission with a tempera- 
ture kToo . Therefore, we use the same spectral distribution 
to describe the emission from the neutron star's surface. 
The angular dependence of the emission is described by a 
Hopf function (Chandrasekhar 1960), which is appropriate 
for the scattering-dominated atmosphere present during a 
burst (e.g., Madej 1991). Photons are propagated to the 
observer through a Schwarzschild metric about a compact 
object with compactness p = R/2M. Corrections to the 
metric due to the unknown density structure of the neu- 
tron star should be of order a few percent (e.g., Braje et 
al. 2000), which is of the same magnitude as the measure- 
ment uncertainties in the observed waveforms. We assume 
that the neutron star has a mass of 1.4 M Q , and denote 
its rotation frequency with Q. 

We consider one single and two antipodal circular bright 
regions with a constant temperature, with the rest of the 
surface taken to emit at a lower temperature. The circu- 
lar hot region(s) are taken to have angular radius p and to 
be centered at an angle a from the rotational axis of the 
neutron star. The angle between the line-of-sight of the 
observer and the spin axis of the neutron star is denoted 
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Fig. 3. — Illustration of the energy dependence of the oscillation 
profiles from our theoretical models of a bright region on a rotating 
neutron star. We model a 1.4 Mq neutron star with a radius of 
10 km (p = 2.5) spinning at f2 = 600 Hz. The bright region has a 
size p = 90°, is located at a = 90°, and is viewed from /3 = 90°. 
The response of the RXTE PCA was accounted for in computing 
the model profiles. The Doppler and time delay effects introduced 
by the rotation of the star cause the amplitude of the oscillation in- 
creases by 6% between 1-20 keV, while the phase of the oscillation 
at 20 keV precedes that at 1 keV by 0.055 cycles. 

by (3. In total, our simulations include seven parameters 
that can affect the energy dependence of the oscillations: 
the compactness and the spin frequency of the neutron 
star; the number, size, position, and temperature contrast 
of the hot regions; and the viewing angle of the observer. 

For each set of parameters, we produce light curves in 
40 phase bins and in 64 energy bins logarithmically spaced 
between 0.01-25 keV. The observed spectrum is simply a 
sum of blackbodies whose temperatures are multiplied by 
Doppler factors; therefore, the signal from a bright region 
of arbitrary temperature can be obtained by rescaling from 
a calculation with kT = 1 keV. For each phase, the result- 
ing spectra are folded through a fiducial PCA response 
matrix, which we generated for PCU 2 during December 
1999 (gain epoch 3), in order to obtain predicted light 
curves that can be compared directly to the observations. 
We analyze these light curves in the seven energy chan- 
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nels used in Section 2, using the same Fourier techniques 
to measure the amplitudes of the oscillations, and a least- 
squares fit to measure the phases. The amplitudes and 
phases we report are measured for the Fourier component 
at the spin frequency of the star when one hot spot is mod- 
eled, and at twice the spin frequency when two antipodal 
spots are modeled. 

3.1. Energy Dependence of the Oscillations 

In Figure 3 we show the effects of the rotation of the 
neutron star on the energy-dependent profiles of the oscil- 
lations {top panel), and on the corresponding amplitudes 
[middle panel) and phases [bottom panel) from a single cir- 
cular bright region of temperature kToo = 2.3 keV at infin- 
ity. Here, the bright region covers half the star [p = 90°), 
is centered at the equator [a = 90°), and is viewed along 
the equator (/3 = 90°). The star has a compactness p = 2.5 
and spin CI = 600 Hz. 

Stellar rotation introduces two main effects: the frac- 
tional amplitudes are larger in the higher energy bands, 
and the pulses at higher energies precede those at lower 
energies [top panel of Fig. 3). Both effects occur because 
the spectrum is steepest at high energies (« 20 keV), so 
that a slight Doppler shift in the apparent photon energy 
results in a significant change in the flux received. This 
not only yields a larger fractional amplitude at w 20 keV, 
but also causes the hard pulse to peak as the spot rotates 
toward the observer, before the observer sees the largest 
solid angle from the spot at lower energies. 

We quantify the energy dependence in Figure 3 in the 
same manner as in Section 2, by measuring the fractional 
rms amplitudes [middle panel) and phases [bottom panel) 
of the oscillations in each energy band. We perform a 
linear least-squares fit to the middle 5 energy channels (2- 



20 keV) from our simulations to determine the slopes of 
the amplitudes and phases as a function of energy. For the 
example in Figure 3, the amplitude increases by approxi- 
mately 0.4% keV -1 , while the phases of the higher-energy 
photons lead those at lower-energy by about 0.003 cycles 
keV" 1 . 

In Figure 4 we display how the amplitudes [top panels) 
and phases [bottom panels) of the oscillations depend on 
the spin and compactness of the neutron star. The left 
panels illustrate the dependencies for the case of a single 
bright region, while the right panels illustrate the case of 
two antipodal regions. The bright regions have a constant 
temperature kT^ = 2.3 keV at infinity and a size p = 60°, 
are located at a = 90°, and are viewed along (3 = 90°. 
The rest of the star is assumed to be dark. 

We find that the energy dependence of both the am- 
plitudes and phases are considerably stronger when the 
neutron star is spinning rapidly [CI = 600 Hz; solid line), 
simply because the Doppler shifts are more pronounced. 
However, the dependence on compactness in Figure 4 is 
more complex. At larger compactness, p > 2.5, the magni- 
tudes of the amplitude increase are very similar for one or 
two bright regions, although the phase lags are about 75% 
larger from two spots. At small compactness [p < 1.7), the 
energy dependence of the pulse profiles from one and two 
bright regions behave quite differently, because the bright 
region is strongly lensed on the opposite side of the star 
from the observer (Pechenick et al. 1983). The lensing 
causes the relative phases at which [i) the observer sees 
the largest solid angle from the emitting region and (it) 
the Doppler effects are strongest to differ for the Fourier 
components at the spin frequency and its harmonic. We 
use p = 2.5 throughout the rest of our simulations, which 
corresponds to a 1.4 M Q star with a 10 km radius. 
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Fig. 5. — Effects of changing the viewing angle f3 (x-axis) and the 
position of the bright region a (solid, dotted, and dashed lines) on 
the slopes of the amplitude and phase as a function of energy. A 
single bright region with p = 60° was modeled, on a p = 2.5 neutron 
star with £1 = 300 Hz. 

We also find that as the size p of the bright regions in- 
creases, the slope of the amplitude as a function of energy 
decreases monotonically. This is due to the facts that (i) 
the amplitudes of the oscillations decrease in general, and 
(ii) points on the bright regions have a range of velocities 
with respect to the observer, which damps the magnitude 
of Doppler effects. The phase lags do not vary as a func- 
tion of the size of the bright regions, because the phase 
of the peak of the pulse is not changed by the effective 
background provided by the regions of a larger spot that 
are not Doppler-shiftcd. 

In Figure 5, we illustrate the effects of changing the 
viewing angle (/?) and position (a) of one bright region on 
the energy dependence of the amplitudes and phases. Here 
we take p = 90° and fcToo = 2.3 keV, and assume that the 
rest of the star is dark. The neutron star has p — 2.5 
and = 300 Hz. We find that the slopes of both the 
phases and amplitudes generally decrease in magnitude as 
the bright region or the observer are moved closer to the 
poles, because the relative velocity of the spot and the ob- 
server decreases. The slope of the fractional amplitude as 
a function of energy decreases more rapidly because the 
amplitude of the oscillation decreases, just as occurs when 
the size of the bright region is increased. 

As the temperature of the bright regions increases, 
the slopes of the phase lags and amplitudes flatten very 
slightly, because the peak of the energy spectrum moves 
to higher energies. However, if the rest of the neutron star 



Fig. 6. — Effects of changing the temperature of the rest of the 
neutron star, while keeping the bright region at a constant temper- 
ature of fcToo = 2.3 keV as observed at infinity. We have taken 
p = a = f3 = 90°, Q = 600 Hz, and p = 2.5. Top panel: The am- 
plitude of the oscillation in the full PCA bandpass decreases mono- 
tonically as the temperature of the surface approaches that of the 
bright region. Center panel: The slope of the amplitude as a func- 
tion of energy reaches a maximum when the rest of the star emits 
within the PCA bandpass, such that at low energies the oscillation 
is muted by a phase-independent background flux, while at high 
energies the emission is only observed from the hot region. Bottom 
panel: Varying the temperature contrast only slightly affects the 
slope of the phase as a function of energy. 

also emits within the bandpass of the detector, the energy 
dependence of the oscillations can become quite dramatic, 
as we illustrate in Figure 6. The slope induced by the tem- 
perature contrast is more than five times larger than that 
produced by Doppler motion alone (i.e., when fcT co i = 0). 
At low energies, the cooler regions on the neutron star 
contribute a constant background flux to the oscillations, 
decreasing their fractional amplitude. At high energies, 
where the cool regions of the neutron star contribute little 
flux, the amplitude of an oscillation is similar to that which 
one would see if the neutron star were dark. The largest 
slope attainable increases as the size of the hot region p de- 
creases. On the other hand, the phases of the oscillations 
are weakly dependent on the temperature contrast from 
the neutron star, for the same reasons that the phases are 
not affected greatly by varying the sizes and positions of 
the bright regions. Thus, the phase lags provide the best 
measure of the apparent velocity of the bright region. 
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4. DISCUSSION 
4.1. Amplitude as a Function of Energy 

In all of the burst oscillations, we find that the ob- 
served increase in amplitude as a function of photon en- 
ergy is about 50% larger than is expected from Dopplcr 
effects alone. In the sources with w 600 Hz oscilla- 
tions, 4U 1636-536 and MXB 1659-298, the amplitude 
increases by up to 0.8-0.9 % keV -1 (see Table 3 and 
Fig. 1), while a single bright region on a star spinning 
at 600 Hz can yield a 0.6 % keV -1 increase, and two an- 
tipodal regions on a star with il = 300 Hz can give rise 
to a 0.3% keV -1 increase (Fig. 4). Similarly, the ampli- 
tudes of the « 300 Hz oscillations from 4U 1702-429 and 
4U 1728-34 increase by 0.5-0.6% keV _1 , compared to the 
0.2% keV -1 increase that might be expected from a bright 
region on a neutron star with Q = 300 Hz. These differ- 
ences are significant, as our uncertainties on the measured 
slopes are less than 0.01% keV -1 . 

It is therefore likely that the observed increase in ampli- 
tude as a function of energy occurs because the rest of the 
neutron star also emits in the the lower energy bands of 
the PC A (Fig. 6), thus reducing the flux variations at these 
energies. To estimate the magnitude of the temperature 
contrast that could produce these oscillations, we use the 
constraints on the sizes and locations of bright regions that 
are consistent with the observed oscillation profiles from 
Muno et al. (2002b). Using the lack of harmonic signals, 
Muno et al. (2002b) show that the pulsations can only be 
produced by (i) a single bright region that covers nearly 
half the neutron star (p « 90°), (ii) a single bright region 
that forms within a < 20° from the rotational pole, or 
(Hi) two antipodal bright regions that form within a few 
degrees of the rotational equator a w 90°. 

A single p = 90° bright region (case i) would produce 
oscillations with an rms amplitude of 55% if the rest of 
the star were dark (Fig. 6, top panel). Since the observed 
amplitudes of the oscillations are 3-5 times smaller even in 
the highest energy band, it is necessary to assume that the 
other half of the neutron star also emits, but at a lower 
temperature, in order to reproduce the observed oscilla- 
tions. A temperature contrast of about 0.15 keV between 
the two hemispheres would produce oscillations with an 
amplitude of 5-10% in the full PCA bandpass (Table 1 
in Muno et al. 2002b) and an amplitude increase of 0.7% 
keV -1 (Table 3). In cases (ii) and (Hi), on the other hand, 
the sizes of the spots are not well-constrained, and thus one 
can reproduce the observed amplitudes with cither small 
hot spots and a small temperature contrast, or with large 
bright regions and a larger temperature contrast (compare 
Fig. 6 with Fig. 4 in Muno et al. 2002b). 

Finally, we note that the increases in the ampli- 
tudes of the burst oscillations as a function of en- 
ergy are not a generic feature of all coherent pulsations 
from LMXBs. The amplitudes of the persistent pulsa- 
tions from the accretion-powered millisecond X-ray pul- 
sars SAX J1808.4-3658 and XTE J0929-314 actually de- 
crease in amplitude between 2-10 keV. This suggests that 
a different emission mechanism or geometry may play a 
role in the generation of pulsations in these two classes of 
LMXBs. 



4.2. Phase as a Function of Energy 

As can be seen from Table 4 and Figure 2, the phase lags 
observed from all of the combined pulse profiles are con- 
sistent either with being constant as a function of energy, 
or with increasing slightly such that the hard pulse arrives 
after the soft pulse. Only one set of oscillations, observed 
from 4U 1636—536 in gain epoch 5, is consistent with a 
negative slope, and in this case an energy-independent 
phase is also an acceptable fit. The magnitudes of the 
phase lags are on order 0.5-0.1 cycles between 2-20 keV, 
which translates to a time delay of « 150 [is for spin fre- 
quencies between 300-600 Hz. Thus, the signs of the ob- 
served phase lags are generally opposite the negative phase 
lags that are expected from Doppler effects, although 
the magnitudes of the lags are comparable. In contrast, 
the sign of the phase lags observed from the accretion- 
powered millisecond X-ray pulsars SAX J1808.4— 3658 and 
XTE J0929— 314 are consistent with those that would be 
expected from Doppler effects, as the pulse at 10 keV is 
observed to precede that at 2 keV by 200-700 ps (Cui, 
Morgan, & Titarchuk 1998; Galloway et al. 2002). 

To understand the sign of the lags, we consider whether 
interactions between photons from the neutron star and 
material (i) in the accretion disk or (ii) in a hot corona 
could modify the energy dependence of the phases of the 
oscillations. In the first case, photons that get scattered 
by the accretion disk can acquire a Doppler shift from the 
motion of the orbiting material. However, since the disk is 
most likely to rotate in the same direction as the neutron 
star, the soft pulse would lag behind the hard (Sazonov & 
Sunyaev 2001), which is still inconsistent with the data. 
In the second case, a corona of hot electrons could in- 
duce hard phase lags if they inverse-Compton scatter pho- 
tons from the neutron star to higher energies (e.g., Miller 
1995). Indeed, the hard X-rays (> 10 keV) observed be- 
tween bursts from neutron star LMXBs are thought to 
originate from such a corona (e.g., Barret et al. 2000). It 
is not clear whether such a corona also exists during X-ray 
bursts, so evidence for it must be sought via simultaneous 
modeling of the burst spectra and the oscillation profiles. 
We also note that a hot corona of electrons could reduce 
the harmonic content of the oscillation profiles (e.g., Miller 
2000; Muno et al. 2002b). 

5. CONCLUSIONS 

We have examined the energy-dependence in the ampli- 
tudes and phases of a sample of 51 burst oscillations using 
data from the RXTE PCA (Table 2). We found that the 
fractional rms amplitudes of the oscillations increase as 
a function of energy by 0.25% keV" 1 to 0.9% keV" 1 be- 
tween 3-20 keV (Fig. 1 and Table 3). We then modeled 
the oscillations as flux variations arising from temperature 
patterns on the surfaces of rapidly rotating neutron stars, 
and calculated the Doppler effects on the energy-resolved 
light curves. Comparing the models with the data, we 
found that the observed slope of the energy dependence 
is generally larger than would be expected from Doppler 
effects alone (Fig. 3). However, it can be reproduced by 
assuming that cooler regions of the neutron star also emit 
in the low energy bands of the PCA (Fig. 6). 

We also found that the high-energy pulses generally lag 
behind the low-energy pulses in the observed light curves 
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by 0.002(1) cycles rad -1 to 0.007(3) cycles rad -1 between 
3-20 keV (Fig. 2 and Table 4). In contrast, Doppler ef- 
fects should cause the high-energy pulses to precede the 
soft pulses by 0.003 cycles keV -1 (Fig. 3). Hard lags like 
those observed could conceivably be produced if the signals 
from the surfaces of these neutron stars are re-processed by 
scattering in a hot corona of electrons. However, it is not 
clear whether such a corona exists during the declines of 
the thermonuclear bursts that produce these oscillations. 
Further observational and theoretical study is warranted 
to address this issue. 
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Table 1 

RXTE/PCA Channel-to-Energy Conversion 



Energy (keV) 
during Epoch 



Channel 


3 


4 


5 


5 


2.3 


2.5 


2.5 


13 


5.1 


5.9 


5.7 


18 


6.9 


8.0 


7.8 


28 


10 


12 


12 


36 


13 


16 


15 


53 


20 


23 


22 



Note. — Boundaries for the 
seven energy intervals used in 
this paper, where the first in- 
terval is 0-5, and the last 53- 
255. See text for explanation 
of gain epochs. 
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Table 2 

Summed Profiles of Oscillations 



Source 


Frequency 
(Hz) 


Gain 
Epoch 


No. 
Osc. 


Total 
Counts 


Bkgd. 
Counts 


4U 1636-536 


581 


3 


6 


312814 


46871 






4 


5 


219356 


23767 






5 


4 


317775 


39453 


MXB 1659-298 


567 


4 


3 


27782 


5761 


4U 1702-429 


329 


3 


7 


547076 


47668 






5 


1 


38883 


5236 


4U 1728-34 


363 


3 


12 


552565 


70491 






4 


6 


291955 


69510 


KS 1731-260 


321 


3 


4 


249986 


43025 


Aql X-l 


549 


3 


2 


120894 


5754 






4 


1 


133784 


6736 



Table 3 

Amplitudes of Burst Oscillations as a Function of Energy 



Source 


Gain 
Epoch 


Constant 

xl 


Linear Trend 
Slope xl 
[% keV- 1 ] 


4U 1636-536 


3 


17 


0.448(6) 


0.9 




4 


16 


0.390(6) 


1.1 




5 


108 


0.889(5) 


5.4 


MXB 1659-298 


4 


5.9 


0.777(18) 


0.7 


4U 1702-429 


3 


48 


0.523(3) 


2.2 


4U 1728-34 


3 


11 


0.263(4) 


1.6 




4 


126 


0.582(6) 


2.0 


KS 1731-260 


3 


14 


0.462(6) 


0.9 


Aql X-l 


3 


8.8 


0.447(8) 


0.2 



Note. — 5 channels were used in computing xi> so there 
were 4 degrees of freedom assuming a constant amplitude, and 
3 assuming a linear trend. 



Table 4 

Phases of Burst Oscillations as a Function of Energy 



Source 


Gain 
Epoch 


Constant 

xl 


Linear Trend 
Slope 
[10~ 3 cyclcs keV -1 ] 


*; 


4U 1636-536 


3 


1.0 


1.92(14) 


0.6 




4 


2.0 


3.57(14) 


0.7 




5 


0.8 


-0.40(8) 


0.9 


MXB 1659-298 


4 


2.6 


7.74(25) 


0.3 


4U 1702-429 


3 


2.4 


1.97(10) 


1.6 


4U 1728-34 


3 


4.5 


3.41(13) 


3.5 




4 


3.5 


4.31(11) 


0.8 


KS 1731-260 


3 


0.5 


1.58(18) 


0.2 


Aql X-l 


3 


0.8 


3.04(29) 


0.5 



Note. — 5 channels were used in computing xt> so t nere were 4 
degrees of freedom assuming a constant amplitude, and 3 assuming a 
linear trend. 



